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Abstract. For light curves sampled on an uneven grid of observation times, the shape of the power density spectrum (PDS) 
includes severe distortion effects due to the window function, and simulations of light curves are indispensable to recover 
the true PDS. We present an improved method for comparing light curves generated from a PDS model to the measured 
data and apply it to a 50-day long RXTE observations of NGC 4945, a Seyfert 2 galaxy with well-determined mass from 
megamaser observations. The improvements over previously reported investigations include the adjustment of the PDS model 
normalization for each simulated light curve in order to directly investigate how well the chosen PDS shape describes the 
source data. We furthermore implement a robust goodness-of-fit measure that does not depend on the form of the variable 
used to describe the power in the periodogram. We conclude that a knee-type function (smoothly broken power law) describes 
the data better than a simple power law; the best-fit break frequency is ~ 10~^ Hz. 



INTRODUCTION AND OBSERVATIONS 

X-ray variability is a prevailing feature of active galac- 
tic nuclei (AGN). With long data sets from monitoring 
campaigns by X-ray missions such as Exosat, Ginga and 
especially RXTE, it became possible to study the shape of 
the power density spectrum (PDS) over many decades of 
temporal frequency. Locally, the PDS is well-described 
by a power law, but globally, a suppression of power on 
long time scales has been found in a number of AGN 
|1, 2, 3, 4, 5]. Because of the direct proportionality be- 
tween the mass and the Schwarzschild radius of a black 
hole, the timescale at which this suppression becomes 
important is expected to obey a linear relationship with 
the mass. Evidence for such a correlation has indeed 
started to emerge 1 6, 7] and can be extended over several 
orders of magnitude to galactic black hole candidates in 
their low/hard states, suggesting that the X-ray emission 
mechanism is similar for both classes of objects. 

Our method for analyzing unevenly sampled light 
curves is based on simulation techniques used by a num- 
ber of researchers 1 6, 8, 9, 10], but introduces significant 
changes. We apply it to the Seyfert 2 galaxy NGC 4945, 
which is unique among the AGN for which a break in 
the PDS has been detected in that the mass of the central 
black hole of 1.4 x 10^ is known fairly accurately 
from mapping the H2O maser emission IIII1 . A precise 
knowledge of the break frequency in this source should 
thus enable us to calibrate the relationship between the 




FIGURE 1. RXTE PCA light curve for NGC 4945 in the 8- 
30 keV band obtained in late 2002. Bin size is variable, with 
an average of 320 s. Error bars (not shown) are typically ±0.3 
counts/s. 



break frequency and the mass of the black hole. 

NGC 4945 was observed by RXTE in late 2002 for a 
period spanning a total of 50 days. A pointing of 1,400 
seconds on average was scheduled approximately every 
6 hours, while for 7 days centered within the total ob- 
servation period the source was monitored intensively. 
Data from the Proportional Counter Array (PCA) were 
reduced using standard RXTE PCA analysis tools. The 
average count rate and the time-averaged photon spec- 
trum are consistent with previous observations lll2L[l3ll . 
The nuclear flux below 8 keV is heavily absorbed at the 



source, presumably by the same material that is respon- 
sible for the megamaser activity. The counts from all 
three layers in PCU and 2 corresponding to the nom- 
inal range of energies from 8 to 30 keV were added to 
produce a combined light curve (16 s intrinsic binning), 
resulting in a total of usable data of about 300 ks. The 
light curve is shown in figure[l] 



DATA ANALYSIS 

The algorithm of Timmer & Koenig fl^ is used to 
generate simulated light curves based on a user-chosen 
PDS model that are then sampled to match the obser- 
vation times of the source light curve. Further process- 
ing comprises subtracting the mean, scaling each light 
curve to the intrinsic variance of the source light curve, 
and adding random Gaussian numbers multiplied by the 
error bars of the source light curve. We then calculate 
the Lomb-Scargle periodogram 

IB in 

and re-bin it 

to 5 points per decade. Using the ensemble of (typ- 
ically > 500) simulated power density spectra, we fit 
the distribution of the periodogram power in each fre- 
quency bin with a stretched distribution: Py{x^,b) = 
b~^py (x^/b) where py is the stretched, py the regular x^ 
distribution, v is the number of DOF, and b is the stretch 
factor. Finally, we compare the source periodogram to 
the fitted distributions, using a ;^^-like statistic based on 
the likelihood ratio as the goodness-of-fit estimator 1 17]. 
The probability that an assumed model with the given set 
of parameters can be rejected is then the fraction of sim- 
ulated light curves that have the likelihood ratio x^ lower 
than the source light curve. The process is illustrated as 
a diagram in Fig. |2l but further explanation is outlined 
below. 

Scaling of simulated light curves to variance of 
source data (step D): In previous implementations of 
light curve simulations, the normalization of the PDS 
model was set to an arbitrary value at first. The best fit 
normalization was then determined by multiplying all 
simulated spectra b^ the same factor until the fit statis- 
tic was minimized |6, 9]. Due to the stochastic nature of 
the periodogram, the light curves simulated from a fixed 
normalization show a spread in variances fisl] . and the 
distribution of the power in each frequency bin will re- 
flect this spread. By comparing the source data to a set 
of light curves simulated in this way, we are effectively 
asking the question: How well do the chosen shape and 
normalization of the PDS model describe the variability 
of the source? In reality, however, we are only interested 
in characterizing the shape of the periodogram. By al- 
lowing the normalization to vary from one light curve to 
the next, a slightly different question can be investigated: 
How likely is it that the source light curve was produced 



by a process that has the chosen PDS shape ? By scaling 
each light curve to the intrinsic variance of the source 
light curve (and adding the Poisson noise level — see be- 
low), we ensure that the area under the curve for each 
simulated spectrum is equal to the one for the source. 

Preservation of the Poisson noise level through the 
normalization of the PDS model (step E): The uncer- 
tainties in the source count rate are expected to con- 
tribute a constant level of power in the periodogram (usu- 
ally called the Poisson noise level). Since the raw count 
rate in the RXTE PCA is dominated by the background 
for this source, the errors on the background- subtracted 
light-curve are approximately Gaussian. To mimic the ef- 
fect of measurement errors in the simulated light curves, 
random Gaussian numbers multiplied by the error bars 
of the source light curve are added. This is done after the 
light curves have been scaled to the correct variance to 
ensure that the Poisson noise level in the simulated light- 
curves equals the one expected in the source data. 

Use of robust fitting statistic to find best-fit parame- 
ters of the model PDS (steps H-J): A useful goodness- 
of-fit measure for non-Gaussian distributions is the rejec- 
tion probability, calculated as the fraction of simulated 
light curves that have a value of the test statistic lower 
than the source data. The choice for test statistic is up to 
the investigator, and an obvious choice is the x^ statis- 
tic, which uses the ensemble average and standard devi- 
ation of the power in each frequency bin for comparison 
to the source power. Because of the non-Gaussian na- 
ture of the distribution of power in each frequency bin, 
the applicability of the x^ statistic is however severely 
limited: specifically, even small departures from Gaus- 
sian distributions have a marked effect on the standard 
deviation and thus on the rejection probability. Since the 
distribution of the power is more symmetric when plot- 
ted in the logarithm compared to the linear power, sub- 
stantially different goodness-of-fit values can be obtained 
simply by a change of variable. The rejection probabil- 
ity, in effect, not only measures how well the model (and 
its associated parameters) describes the source data, but 
also the degree of non-Gaussianity. In contrast, by fitting 
the distributions of the power in each frequency bin with 
suitably parameterized functions (step H in the analysis), 
the effect of variable transformation on the fitting statistic 
is removed: If stretched x^ distributions are an adequate 
description of the distribution of linear power, then under 
a change of variable the data can equally well be fit with 
the similarly transformed functions. The probability den- 
sity for observing the source power as well as the value 
at the peak of the distribution is unchanged under these 
transformations. The likelihood ratio is thus an invari- 
ant quantity, i.e. independent of the chosen form of the 
variable used to describe the power in the periodogram. 



Analysis 



Description 



^ step 

^ starting with . . 



■p> generating 
light curves 



rebinning of 
the light curve 



RXTE lightcurve, 

16 sec intrinsic binning, 

background-subtracted 



average newbin size 

320 sec, variable between 

256 - 384 sec 



Simulated data 

PDS model, values for 
model parameters, 
normalisation arbitrary 



use Timmer & Koenig 
[1995) algorithm, 16 sec 
binning, contiguous, 
length equal to source 
light curve 



sampled and rebinned in 
same manner as source 
light curve 



Y\ variance 
adjustment 



E 



measurement 
errors 



calculation of Lomb- 
Scargle periodogram^*^ 



error bar variance 
subtracted from total 
variance to determine 
intrinsic variance' 



light curve mean- 
subtracted and scaled to 
intrinsic variance of 
source light curve 



add random Gaussian 
numbers multiplied by 
error bars of source 
lightcurve 



evaluated on even grid of frequencies up to 
1 / (2 X 320 sec) 



G 



rebinning of the 
periodogram 



TJ fitting the distribution 
of power 



T likelihood ratio as 

goodness-of-fit 



rebinned to 5 points per decade in the logarithm 
of the frequency 



' use ensemble of 
simulated spectra, fit 
distribution in each 
frequency bin with x^- 
distributions with two 
adjustable parameters: 
no. of dof and stretch 
factor* 



repeat for source spectrum and all simulated spectra, 
read off probability density for observing given power 
[using above parameterisation), divide by maximum of 
distribution, sum up logarithm of (-2 x ratio) for all 
frequency bins, call this sum the likelihood ratio 



fraction of simulated lightcurves that have likelihood 
tj rejection probability ratio %^ lower than for the source lightcurve is probability 

that assumed model with current values for parameters 
can be rejected 



FIGURE 2. The diagram of the analysis steps used in the process of PDS characterization. 
Footnotes to the diagram: 

(1) error bar variance is calculated as the mean square of the source light curve error bars 

(2) see Lomb 1 15] and Scargle 1 16] 

(3) Since the rebinned light curve includes contiguous segments sampled at 320 sec resolution in the intensive section of the 
observation, information about the PDS shape can be extracted down to that timescale. 

(4) We are only concerned here with finding a suitable parameterization of these distributions, distributions with a stretch factor 
are a natural choice since the power in the periodogram is drawn from distributions scaled by the power of the underlying PDS. 




FIGURE 3. Rejection probability for the unbroken power FIGURE 4. Rejection probability for the knee model, show- 
law model, The best fit slope is 1.3, with a rejection probability ing contours for 50, 68, and 90%. The best fit is obtained with 
of 70% indicating a poor fit. a break frequency = 10"^ Hz and a slope a = 1.5. The 

corresponding rejection probability is 40%. 



RESULTS 

We investigated two models: the simple (unbroken) 
power law (power law slope as only parameter), and 
the knee model of Uttley, McHardy & Papadakis |9] 
(p{f) = ^[1 + {f / fb)^]~^^^'^ parameters: high frequency 
power law slope a, break frequency ft). Figures 13 and|4] 
show the rejection probabilities obtained. The unbroken 
power law results in a minimum rejection probability of 
70% at a slope of 1.3, indicating a rather poor fit. The 
rejection probability is reduced to 40% by introducing 
a 10~^ Hz break at a high frequency slope of 1.5. It is 
clear that the exact location of the break is dependent on 
the particular PDS model used to characterize the source 
variability. 

Uncertainty on rejection probability: To quantify 
the effect of simulating a finite number of light curves 
from which the goodness-of-fit statistic is estimated, we 
simulated more than 800,000 light curves using the un- 
broken power law model with a slope of 1 .3. The process 
of fitting the distributions of power through calculating 
the rejection probability (steps H through J in Figure |2|l 
was repeated for sets of light curves of increasing size. 
The mean of the distribution of rejection probabilities is 
consistent across the different sizes, while the standard 
deviation decreases from about 3% to 2.5% by going 
from 500 light curves to over 16,000. Assuming that the 
spread of rejection probabilities is only a weak function 
of the model and its associated parameters, we quote a 
characteristic uncertainty on our fitting statistic of about 
±3%, based on 500 simulated light curves. 

Simple power law vs. knee model: The difference 
in rejection probability between the unbroken power law 
and the knee model is much larger than the above uncer- 
tainty. The flattening of the power towards low frequen- 
cies in NGC 4945 has thus been detected at a statistically 
significant level. The length of our observation is most 



likely not sufficient to distinguish between models that 
use different functional forms for the break in the PDS. 
In particular, the differences in rejection probabilities be- 
tween the knee model and the broken power law model 
used in Uttley, McHardy & Papadakis 1 9] is not expected 
to be statistically significant for the currently available 
light curve of NGC 4945. 
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